Spin relaxation and antisymmetric exchange in n-doped III-V semiconductors 
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Recently K. Kavokin [Phys. Rev. B 64, 075305 (2001)] suggested that the Dzyaloshinskii-Moriya 
interaction between localized electrons governs slow spin relaxation in n-doped GaAs in the regime 
close to the metal-insulator transition. We derive the correct spin Hamiltonian and apply it to the 
determination of spin dephasing time using the method of moments expansion. We argue that the 
proposed mechanism is insufficient to explain the observed values of the spin relaxation time. 
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Observation of exceedingly long spin relaxation times 
for electrons optically pumped into the conduction band 
of n-doped GaAs Q prompted suggestions for the use of 
polarized spins in quantum computing . No consensus 
has been reached, however, on the mechanisms govern- 
ing the spin relaxation processes. Low-temperature spin 
dephasing time r s in Ref. jl| behaves non-monotonically 
as a function of donor concentration n d , reaching values 
above 100 ns at rip — 10 16 cm -3 . This concentration 
still belongs to the nonmetallic regime. 

In an interesting paper || the author introduced a new 
dephasing mechanism acting between the spins of elec- 
trons localized on randomly distributed shallow donors. 
According to Ref. jj), dephasing would occur each time 
the spins of the two electrons on neighboring centers ro- 
tate against each other in an effective spin-orbital field 
in the process of tunneling between the two centers dur- 
ing the spin flip-flop act caused by the exchange inter- 
action. This notion and the way it was treated in Ref. 
S in terms of the exchange narrowing does not seem 
to be fully self-consistent. The point is that "random 
effective fields" originating in pjfrom the antisymmet- 
ric Dzyaloshinskii-Moriya (DM) [|, || interaction and the 
dynamical averaging due to the isotropic exchange in fact 
both arise from the interaction of the spins of the very 
same pair and hence are not independent. 

As a more consistent approach to the problem we tried 
its formulation in terms of the EPR line shape for the 
spin Hamiltonian which for semiconductors without an 
inversion center is a sum of the conventional isotropic 
and the antisymmetric DM exchange terms. The latter 
is responsible for non-conservation of the total spin of the 
system. 

The antisymmetric interaction 



Hum = d(R)(§i x s 2 ) 



(1) 



in semiconductors lacking the inversion symmetry (III- 
V and II- VI compounds) is of interest in itself, but has 
not been properly calculated yet. In Ref. || the asymp- 
totic behavior of d(R) at large distances between centers 
R ^> ob, where ae is the effective Bohr radius, for shal- 
low (hydrogen-like) donor centers was calculated in the 
Heitler-London approximation. But it has been shown 



H [T(| that the Heitler-London theory does not account 
for the Coulomb correlations inside the "tunneling corri- 
dor" connecting the two centers which dramatically af- 
fects the prefactor of the isotropic exchange. 

Below we derive the asymptotic form of d(R) at R ^> 



ae- We have slightly modified the approach of Ref. j£ , 10 
to include the initial spin degeneracy of the problem and 
its lifting by interactions in the presence of a weak spin- 
orbital interaction 



h(p)s. 



(2) 



Our theoretical result for the -R-dependence of the an- 
tisymmetric exchange d(R) differs significantly from Ref. 
H 0, H . In Ref. [Q the authors concluded that their ex- 
perimental results for the spin relaxation time t s together 
with the data || agree well with the mechanism || in 
the concentration range nry — 2 x 10 15 -f 2 x 10 16 cm -3 . 
Upon substituting into the same formulae of @], our 
corrected DM term increases the estimate for t s by a fac- 
tor of 4. Although this does not rule out completely the 
relaxation mechanism || because of the phenomenologi- 
cal character of the analysis of the experimental data in 
, the possibility that the antisymmetric exchange alone 
accounts for the values of r s observed in |l|, |[ fjj in this 
doping interval seems rather unlikely. 

The essence of the method is as follows. The envelope 
of the wavefunction of an electron localized on a donor 
impurity in the effective mass theory satisfies Schrodingcr 
equation with the Hamiltonian 



H = p 2 /2m e + H s 



(3) 



where m e is the effective electron mass, n is the permit- 
tivity. (For GaAs m e rj 0.072too and k « 12.5). 

The spin-orbital Hamiltonian (^|) in zinc-blende semi- 
conductors has the form [1TT1 FHl 



Mp) 



m e y / 2m e E, 



--Px(p 2 y 



vl). 



(4) 



h y (jp) and h z (p) are obtained by permutation in (|j). 
Here E g is the band gap and a so is a phenomenologi- 
cal parameter. For GaAs E g ss 1.43 eV and a so » 0.07. 
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Without the spin-orbital interaction the ground state 
had the energy ^E = — Ry = —m e e A /2k 2 . (Eq rs —12.6 
meV for GaAs). The corresponding eigenfunction of (|^) 
were <p^(r) = cp (r)5™, where = 5?, denotes the 
spinors for up and down spin projection, 



1 



-r/oB 



(5) 



and Ob = K/m e e 2 is the Bohr radius (ae ~ 92 for GaAs). 
In the presence of a spin-orbital interaction the spin pro- 
jection no longer commutes with (||), but the two-fold 
(Kramers) degeneracy of the ground state remains. For 
small spin-orbital interaction (|j) eigenfunctions of (||) are 



the function tp^ should convert to (p"(ri&)(pP(r2B), i.e. 
ycta ,pp _^ fiaa fi@0 _ rp^ e p rmc jp a i terms in spin-orbital 
interaction are included into definition (^|) , and the equa- 
tion for x aa ,/3 ' 9 turns out to coincide with the one in the 
absence of the spin-orbit 0] 



dr, -d. 



1 



1 



i2 — 2Ci R + R ria 



x = o. 



(10) 

Here we used that dip ~ <p/clb and hence one can neglect 
all but the first derivatives of the functions x which varies 
on distances of order R 3> ob ■ 

Therefore \ aa ,/3 ^ = xoS aa S 13 ^ , where on the median 
5-dimensional hyperplane £1 = (2 = C @ 



<(r) 



^ (r)(e imeah ^ r<T / 2 °B) c 



(6) 



where f = r/r. The phase factor turns into unity when 
r — > 0. Thus the Kramers index /1 takes on the meaning 
of a spin projection of an electron on the center. 

The two-electron Hamiltonian is obtained by com- 
bining two single-electron Hamiltonians (0) and the 
Coulomb interaction of electrons between themselves and 
with the alien ions: 



Utot = (Pi + P 2 ) / 2 ™e + Wsol + Ko2 

f! f J_ I JL J_ J_ 

« \ r 12 -R riA riB ^2A 



(7) 
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Here riA,B and ?"2a,b denote the distances between the 
electrons and the two donor centers A and B. 

For isolated donors, i.e. on large distances R ^> clb 
between two centers a two-electron wavefunction is 
merely a product cp" (ri a) ^(t2b) of two one-electron 
wavefunctions (||). For large but finite R these func- 
tions have to be corrected to account for potentials of 
the Coulomb interactions of electrons between themselves 
and with the alien ions (0). Following one starts 

with the Schrodinger equation for the corrected functions 



njaia',/30' I a' j3' 
rl tot YthV 



(8) 



Let C denote the direction from the donor A to the 
donor B, positioned at £ = ^fR/2 respectively. It turns 
out H that for the calculation of exchange interaction it 
is enough to know ip^ on the median hyperplane £1 = C2 
far from the centers. There one may seek for the func- 
tions ip^ v approximately in the form 

rj(ri,r 2 ) = x Qa '^Vi,r 2 )^'(r 1A )^'(r 2B ), (9) 

where x Qa ,/3 ' 9 varies on the scale of the order of R. 

Substituting into (@) yields a differential equation 
on x aa ■ The boundary condition is determined from 
the condition that when either i^a — > or r2B — > 



2R 



2pi 



i? + 2|C|V^-2|CI 



e 2R . 



(11) 



Here p\2 = \pi — p 2 l) and Pi 2 are the transverse radii- 
vectors of the electrons. 

The two-electron eigenfunction i/> a/3 (ri, r 2 ) of 7i to t ad- 
justed for the exchange is then sought in the form 



4, 



a/3 



r %b ap + 



'run J 



(12) 



where and cja, are unknown coefficients. VC/t stands 
for the interchanged states of the two electrons, when the 
first electron is now almost localized on atom B in the 
state v and the second — on atom A in the state /i 



C^(ri,r 2 )^ Q ( r 2,ri) 



(13) 



Because the total fermionic wavefunction should be 
antisymmetric in particle permutation ip al3 (ri, r^) = 
— ^^ Q (r 2 ,ri), the coefficients effi = — c^ v . 

We can now proceed to the derivation of the total ex- 
change Hamiltonian. Multiplying the Schrodinger equa- 
tion Eip = Tiip by tp* (for brevity we omit spinor struc- 
ture), and the Schrodinger equation (||) H*ipfiu* = EqiJj^ 
by ip, subtracting term by term and integrating over the 
region £1 < C2, we find 

(E - E ) Cflv = f fe,#V " W*^) rf 6 ri, 2 . (14) 

The kinetic part of Ti. can be readily reduced to a sur- 
face integral || 10 



(E - E )c^ = \ j (Wp, - ^diP) dS (15) 

over the median plane. Here dS = — £ 2 )d( t d' l p 1 2 
is its surface element and d = (81,82) is the six- 
dimensional gradient. 

We then substitute (|l^) into ([l5|) neglecting derivatives 
of Xo compared to those of ip^ : 



(16) 
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After simple calculations Eq. (15) results in 



{E - E )c^ 



-2c 



(P) 



Jiv/2s (P) 



= J(R)(e-^/ 2 )^(e 

with the angle 

7 = m e hh(C)R/a B 
The exchange integral equals ||, |l0) 



(18) 



J(R) = -2 1 x%cpQ(r 1A )cpQ(r 2B )cp Q (r2A)¥>o(riB)d(d 4 p 12 

Cl=C2 



-0.818 x 2Ry (i?/a B ) 5/2 ^" 2fl / aB 



Using 

we finally arrive at the equation for the coefficients c^, 
(E — Eq)c = 7i e xC, with the Hamiltonian 

-«"V 2 (*i6- 2 + 1) (21) 

CTl i? (7)<T2 + 1 

1 -7(d-l X <X 2 )] , 



« — ij[oi<T 2 

where the Pauli matrices <xi and <x 2 are supposed to act 
in the space of the Kramers indices, and R (7) is the 
three-dimensional rotation matrix on the angle 7 around 
the axis 7. 

The rotation angle 7 ([II]) gives the strength of the DM 
term relative to the isotropic exchange. It is perpendic- 
ular to axis C since h(£)£ = 0. Taking the square mean 
of h(£) over the directions of C we find that for GaAs 



2a s 



R 



7=(7 2 ) 1 . /2 = 

t y/35E g /Ry a B 



0.00157—. (22) 



This 7 is shown in Fig. [I] against the i?-dependence of 7 
from Ref. fft Eq. (7)] for comparison. In the region of 
interest (R ~ 3 -7- 7a B ) our result for 7 is about one half 
that of Ref. @, Eq. (7)]. 

We turn now to the question of how the antisymmetric 
exchange would manifest itself as a spin relaxation mech- 
anism. It was suggested in Ref. || f?J to write down the 
corresponding relaxation time T sa in terms of a precession 
mechanism |13|: 



2 ' c ' ' 



(23) 



where now, unlike ]l2] ] , l/r c is a rate of the flip-flop tran- 
sition caused by isotropic exchange. Inconsistency of the 
physical picture of Ref. || in the use of the dynamical 
averaging in Eq. ( pSj ) was mentioned above. Another 
objectionable point is to write r c in the form |Q 



as 
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FIG. 1: Mean square rotation angles 7 of the asymmetric in- 
teraction in GaAs as functions of the inter-center separations 
as calculated in the present paper (thick line) and according 
to the interpolation Eq. (7) in Ref. Q (thin curve). 



where R c ~ , £ and (3 being the fitting parameters 

of the order of unity. The choice for the latter one, (3 w 
0.65, interpolates the maxima in the Poisson distribution 
for the nearest and next nearest neighbors. Such a form 
is assumed in Ref. |?J to phenomenologically account for 
large average distances between donors. But since the 
maximum in the Poisson distribution is rather broad, it 
appears erroneous to merely substitute R c for R in the 
exponential dependence J(R). 

At no ~ 10 16 cm" 3 the distance R c « 3.3ob, and 
from Fig. |l| for the angles 7 |7| Eq. (7)] one obtains 
7(i? c ) « 0.009, while the corrected value in Eq. (j||) 
7 rj 0.005. The additional factor « 4 in Eq. ( p3[ ) imme- 
diately erodes the claimed quantitative agreement with 
experimental data in Ref. jl], 

With the applicability of (^) being questioned and the 
experimental value of t c being not determined in un- 
ambiguously, we try to explore spin relaxation due to DM 
terms as a question of the EPR line shape for a dilute sys- 
tem of paramagnetic spins interacting via Hamiltonian 
( pl| ) . This appears a proper approach if all other sources 
of inhomogeneous line broadening are neglected and the 
limit is taken of the Larmor frequency wl = — > 0. 
In Ref. p 



h/£J{Rc 



such a problem dealing with DM inter- 
action has already been addressed for Mn-based II- VI- 
compound diluted magnetic semiconductors in the frame- 



(24) work of high temperature moment expansion of the linear 
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response function 



e- iut ([M x (t),M x (0)})dt. (25) 



Here square brackets denote quantum commutator of 
the magnetization operator M = g/i B J2 n s ™ (summa- 
tion runs over all donors). Angle brackets denote both 
thcrmodynamical average and averaging over disorder. 
The spin relaxation time is then given by p3] 



2I 4 /wli, 



(26) 



where the frequency moments of the response function 
I n = ^ oo Lo n - 1 x"{oj)duj/TT X {0) obtained in Ref. g| as 
a high-temperature expansion, are 

h = -Tr([7i DM ,Mj 2 )/Tr(M2) + 0(T- 1 ), (27) 
h = Triin^iH^AQ^/TviM^ + OiT- 1 )^) 

Calculating the commutators we find 

h = ±n D J J 2 (r)[ 7 2 (f)- 7 2(f)]d 3 r + 0(T- 1 ),(29) 

h « 32ns j J 4 (r) [d^W-jKr^cPr + OiT-tyO) 

(For an order of magnitude estimate below we neglected 
in (30) the terms with double integration under the pre- 
text that they are of the "second order" in rijj. Averaging 
over f we use {ji(r)jj(r))? 



The high-temperature expansion (29), (pO]), where one 



would integrate over all randomly distributed centers, ap- 
plies to weakly interacting spins only. Eqs. (|2^), (|30| ) 
need to be corrected to take into account that if two spins 
are so close that their exchange interaction J{R) ^> T, 
they would form a singlet state and drop out of the 
system's thermodynamics. A proper calculation of the 
temperature effects in the EPR line due to clusterization 
should proceed in the framework of a scaling theory p3 . 
The problem is difficult and remains unsolved. We take 
the effects into account qualitatively introducing a cut- 
off in the integrals ( p9| ) , (|3(]) at short distances ar found 
from the condition J(ar) = T. (At T = 2K this radius 
aT ~ 3.7a B , at T = 4.2 K cit ~ 3.1a B ). 

Calculation of (29), ( |30| ) yields the answer in terms of 
incomplete gamma functions T(n,x): 
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'2r(13,8a T /a B ) 



7mD«B0.818Ry7 2 V 7rr 3 (8, 4a T /a B ) ' 



(31) 



i 16 cm" 3 at T = 2K we find r sa 



For GaAs at n D = 10 1 
290 ns, and at T = 4.2K we find r sa « 170 ns as compared 
to t s « 80 ns found experimentally. (Using expression 
( |23| ) with corrected value for 7 from Fig. [l] would give 
the temperature-independent value of r so ~ 330 ns at 
n D = 10 16 cm" 3 ) 



The method of taking a spin-orbital interaction into ac- 
count, developed in the present paper to shallow donor 
centers in bulk zinc-blende semiconductors, may be eas- 
ily generalized. In a low-dimensional environment, e.g., 
in quantum dots, spin-orbital interaction is described by 
a linear dependence on p: hi(p) = hikPi, where hik is a 
tensor. The spin structure of an exchange Hamiltonian 
( |2~l| ) acting on the coefficients of expansion (|l^) again de- 
couples from the coordinate dependence J(R). Although 
explicit expression for J(R), of course, now depends on 
the potential of the quantum dots, the form of the "ro- 
tated" exchange Hamiltonian (|2l]) with the rotation an- 
gle 7^ = m e hhik(kR/aB instead of ( |l8| ) remains. For a 
specific choice of the linear in p spin-orbital interaction 
this result may be obtained by an exact unitary transfor- 
mation as has been recently shown in [fl6|| . 

In the end, we comment on the role of the spin- 
orbital interaction of electrons with the electric field of 
alien ions. For a homogeneous field it has the form 
]l5[ a£/iB<?s(p x ~E)2c/E g , the coefficient oie ~ 1 de- 
pends on the band structure. This term is proportional 
to ~ a£;(Ry 2 / E g )(m e /mo) , and is much smaller than 
the term (Q) considered above that is proportional to 
-a so Ry(Ry/^ 9 ) 1 /2. 

To summarize, we derived the asymptotically cor- 
rect form of the spin Hamiltonian for two hydrogen-like 
donors in the n-doped GaAs. We applied the EPR line 
shape formalism to analyze spin dephasing times T sa due 
to the antisymmetric DM exchange. The rough esti- 
mate for r sa (neglecting the detailed low-temperature 
spin clusterization which would further increase r so ) 
yields values exceeding those observed experimentally. 

— 1/3 

Although the concentration range n D ~ a B is very dif- 
ficult for theory, our results suggest that the anisotropic 
exchange between the localization centers is not a pre- 
vailing mechanism of the spin relaxation in this concen- 
tration range. 
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